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We have calculated the energy per site for the ground state of antiferromagnetic quan- 
tum spin chain with variable range exchange h(j — k) oc sinh 2 a sinh -2 a(j — k) in the 
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framework of the asymptotic Bethe ansatz. By expanding it in powers of e 2ct , we have 
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exchange obtained earlier in the atomic limit of the Hubbard chain. 
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At present, there are two main approaches in investigations of spin correlations in the 
ground state of the quantum Heisenberg antiferromagnet. The first one is connected with 
analysis of the representations of quantum groups [1] and gives some universal prescription 
of calculation of various correlations via multiple contour integrals. Unfortunately, till now 
there are no possibilities to reduce these multiple integrals to known transcendents. The 
second approach consists in finding correlators in the form of Fredholm determinants 
[2]. However, within this approach there is still no way to analytic evaluation of these 
Fredholm determinants. Asymptotics of spin correlations at large distances between spins 
has been obtained in the framework of conformal field theory [3], but this theory cannot 
give the values of correlators at finite distances. 

The first explicit calculation of nearest-neighbor spin correlator < OjOj+i > in the 
thermodynamic limit has been made by Hulthen [4] with the use of the Bethe ansatz [5]. 
Since this correlator coincides with the energy per site in the vacuum state of the model, 
the problem is equivalent to calculation of this energy for finite number of lattice sites N 
and taking thermodynamic limit N — > oo. The second-neighbor correlator < > 
has been calculated by Takahashi [6] with the use of perturbation theory in the atomic 
limit of the half-filled Hubbard model. This result has never been confirmed by other 
scheme of calculations. 

In this letter, we perform the explicit calculation of < OjOj+i > by using the way which 
differs substantially from the way of ref. [6]. Namely, we start from the consideration of 
the integrable spin chain with variable range exchange and the Hamiltonian [7] 

^ _ 1 y-v sinh 2 a a fit - 1 
2 - +k sinh 2 a(j - k) 2 

In the nearest-neighbor limit a — > oo one can decompose (1) into a series 

n = \ EOWi - 1) + ^- 2a £(^+2 - 1) + (e- 2 *). (2) 
j j 

Hence one can write the ground-state energy per site as 

e = \ < a 3 a J+1 - 1 > +h~ 2a < a 3 a J+2 - 1 > +o(e~ 2a ) : (3) 

where < > means average on the vacuum state of the Hamiltonian (2). Fortunately, in 

2 



the first order approximation (3) one can replace this state to the vacuum state of nonper- 
turbed Hamiltonian with the interaction of nearest-neighbor spins, TCq = \ X^^j^j+i - !)• 
Let us now calculate the ground-state energy per site (3) of the model with the use 
of the asymptotic Bethe ansatz (ABA) [8]. The wave functions of the states with M 
down spins can be calculated exactly [7]. Their asymptotic expression in the region 
rt\ Cn 2 ... <C Um has the form similar to Bethe one, 

/ m \ fi m \ 

?/>(ni, ..n M ) oc ex P *E k P*n a exp - x(k Pa , k P/3 ) , (4) 

P67TM V a=l ) \ l a<f3 J 

where the first sum is taken over all permutations from the group ir M , {k a } is the set of 
pseudomomenta and x(k a , hp) is the two-magnon phase shift defined by the relations 

cot^^- = V (K)- V (k l> ), (5) 



2nia \2aJ 2ia \2a / 
where Ci is the zeta Weierstrass function defined on the torus T 1 = C/(Z + ^ : Z). To 

consider the chains of finite length N in the thermodynamic limit iV — > oo, we adopt the 
main hypothesis of ABA, i.e. imposing periodic boundary conditions on the asymptotic 
form of the wave function (4) (it is worth noting here that the exact solution of the problem 
of finite chains in the model should be based on the treatment of the Hamiltonian with 
the exchange given by the Weierstrass p function instead of (1), and is not available till 
now except the cases M = 2,3). Taking ip(n,2, --Um, ni + N) = ip(ni, ..tim) and calculating 
the both sides with the use of the formula (4) results in the ABA equations 

/ M \ 

exp(ik a N) = exp U ^ x(K, M , a = l,..M. (7) 

V ^ J 

The energy of corresponding configuration is given by 

sinh 2 a 
sinh 2 an 



a=l n^O 



We are interested in the antiferromagnetic vacuum of the model and should take iV even, 
M = N/2. Taking logarithms of both sides of (7) and choosing the proper branches [9], 
one arrives at 

= ~~2^ ~VN^ arctan [v 9 ( A; «) - <P( k p)]- ( 9 ) 



We adopt usual hypothesis about the distribution of the (half) integers {Q} for antifer- 
romagnetic vacuum state. Namely, it will be assumed that these numbers form uniform 
string from —Q ma x to Q ma x, Qmax = N/4 — 1/2 without holes. We introduce the rapidity 
variable A by the relation A = (p(k) and the function //(A) via the relation ir — k = //(A). 
The ABA equations (9) now can be written as 

Q a /N = Z(X a ), (10) 

where 

Z(X) = (27r)"V(A) - 4? E arctan(A - A^). 

Following Hulthen [4], let us go to continuous variable x = Q a /N in the limit N — > oo 
and introduce the root density cr N (X) by the relation a N {X) = dx/dX. Differentiating 
both sides of (10) with respect to A, one arrives at the Hulthen-like equation in the limit 
N ^ oo 

^(A) = (2tt)-V(A) - ^ A(X - A')(7 00 (A , )dA / , (11) 

•J — oo 

where A(X) = + A 2 )] -1 . The energy per site can be written from (8) as 

/oo 
e(A;(A)) ( T 0O (A)dA, (12) 
-oo 

where 

l9 JS, cosnk(X) — 1 . , 

e(k(X)) = 2sinh 2 a V -V • (13) 

n=1 sinh an 

The solution to (11) can be found via Fourier transform, 

/oo p^Pflr) f°° 

Substituting it into (12) yields 

/oo roc pip^ roo 

dXe(k(X)) / dp— — - / //(tK^t. 
-oo J — oo 1 + 6 <p< J —oo 

Choosing variables as A = (p(k), fj/(r)dT = —dk' and changing the order of integration (it 
is allowed since the integral over r vanishes sufficiently fast as \p\ — > oo), one arrives at 
the main formula 

/oo Aj\ f2ir r2n 

.ooT+e-W Jo dk£ (W(k)e^j Q dk'e-^ k \ (14) 



where the functions e(k) and tp(k) are determined by (13) and (6). The integrals in (14) 
cannot be calculated analytically. However, one can see from (6) and (13) that for large 
a 

ip(k) = - cot - + 2e~ 2a sin k + o(e" 2a ), (15) 
e(k) = 2(cos k - 1) + 2e" 2a (cos 2k - 1) + o(e~ 2a ). (16) 

Substituting (15) and (16) into (14), one can calculate inner integrals up to the order of 

e~ 2a , 

dk'e~ ip ^ = 2n(l-p 2 e- 2a )e-\ p,2 \ + o(e~ 2a ), 

J o 

/ dke{k)<f'{k)e iMk) = tt(2 + e- 2a (8 - p 2 )) e - |p/21 + o{ e - 2a ). 
Jo 

Now, taking into account the formula J Q °° = 3/2£(3), where ( is the Riemann zeta 
function, the two leading terms of are calculated. The comparison of (3) and (14) 
then gives 

< BjBj+2 >= 1 - 16m2 + 9C(3), (17) 

which coincides exactly with the expression given by Takahashi [6]. 

To conclude, we would like to point out that our approach to the calculation of the 
correlator is based on the hypothesis of the validity of the ABA and uniform distribution 
of the numbers Q a for the antiferromagnetic vacuum of the model (1). The approach of 
ref . [6] refers to the perturbation theory in the atomic limit of half-filled Hubbard model 
and the integral representation of its ground-state energy in the thermodynamic limit 
which is derived under the assumption on uniform distribution of integers in the Bethe 
ansatz equations for the model. Hence both results are not absolutely rigorous from math- 
ematical point of view, and their coincidence supports the validity of both approaches. 
Unfortunately, neither our method nor the method of Takahashi [6] can be used for the 
calculations of higher correlators. It seems that these calculations will be possible if one 
will find the method of reduction of multiple integrals in general formulae of [4] . Till now, 
there is no recipe of doing this even for the second- neighbor correlator (17). 
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